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Abstract 

Shear and bulk viscosities of deconfined gluonic matter are investigated within an effective kinetic 
theory by describing the strongly interacting medium phenomenologically in terms of quasiparticle 
excitations with medium-dependent self-energies. We show that the resulting transport coefficients 
reproduce the parametric dependencies on temperature and coupling obtained in perturbative QCD 
at large temperatures and small running coupling. The extrapolation into the non-perturbative 
regime results in a decreasing specific shear viscosity with decreasing temperature, exhibiting a 
minimum in the vicinity of the deconfinement transition, while the specific bulk viscosity is sizeable 
in this region falling off rapidly with increasing temperature. The temperature dependence of 
specific shear and bulk viscosities found within this quasiparticle description of the pure gluon 
plasma is in agreement with available lattice QCD results. 

PACS numbers: 12.38.Mh;25.75-q;52.25.Pi 
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I. INTRODUCTION 



Transport properties of strongly interacting matter, encoded in coefficients such as shear 
and bulk viscosities, which describe the hydrodynamic response of the system to energy 
and momentum density fluctuations, are of particular importance for understanding the 
nature of QCD matter. Their firm knowledge, besides other characteristics like the QCD 
equation of state (EoS), has great impact on a variety of physical phenomena in cosmology, 
astrophysics or nuclear physics. In particular, the dynamical description of the medium 
created in relativistic heavy-ion collision experiments requires a precise determination of 
the temperature T and density dependencies of these transport coefficients. Even though 
tremendous efforts have been invested both from theoretical as well as experimental sides, 
this remains still a challenging task, especially in the regime of experimentally accessible 
temperatures and densities. 

It was shown that the observed collective flow behaviour of the matter created at the 
Relativistic Heavy Ion Collider (RHIC), in particular the measured elliptic flow jj-yj, could 
quantitatively be described by means of ideal hydrodynamics |5-10] for not too large trans- 
verse momenta of the particles. Although the results of these calculations depend crucially 
on the employed freeze-out prescription as well as the assumed initial conditions 



11 



18|, the 



success of ideal hydrod yna mics suggested at most small dissipative effects in the medium 
produced at RHIC 1 9M2 11] . This observation led to the interpretation that the produced 



quark-gluonp 



acting gas 



22 



asma (QGP) behaves as a strongly coupled fluid rather than a weakly inter- 



24|. Detailed ca 



dynamics were performed 



13 



culations within the framework of relativistic dissipative fluid 



25N38| , confirming indeed the smallness of dissipative effects. 



39 



40[ are solvec . 



41 



43|. 



Therein, the equations of (for many EoS) causal Israel- Stewart theory 
Recently, however, corrections to this second-order in gradients theory were reported 

On the other hand, some transport coefficients, in particular the shear viscosity rj, cannot 
be arbitrarily small. In fact, the quantum mechanical uncertainty relation sets a fundamental 
lower bound on rj for any physical system 44J. Based on unitarity arguments 45M47I] . the 
dimensionless specific shear viscosity (i. e. the ratio of shear viscosity to entropy density 
s) was recently conjectured to be bound from below by {ji/s)kss > h/(ATt) (Kovtun-Son- 
Starinets bound). The comparison of causal viscous hydrodynamic calculations with data 
on low- momentum elliptic flow spectra measured at RHIC showed that rj/s should be close 
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to this lower bound 
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21 
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33 



48| . Nevertheless, some other results seem to indicate 
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361 ] . Moreover, there 



that the KSS bound might be violated for the QGP at RHIC 
are other physical systems 49145 21]. where this is indeed the case. 

The specific shear viscosity ij/sasa function of temperature is known to exhibit a min- 



imum in the vicinity of a phase transition for a variety of liquids and gases 47], |53j as well 



as for ultraco 
conjecture [20 



d fermionic systems close to the unitarity limit 54]. This gave rise to the 



53| that rj/s for QCD matter exhibits a similar behaviour in the vicinity of 



the deconfinement phase transformation. 

In contrast, the bulk viscosity ( is exactly zero for conformally invariant systems. This is 
due to the vanishing trace of the energy-momentum tensor in conformal theories. In QCD, 
as quantum fluctuations break the scale invariance, the bulk viscosity is non-zero at finite 
T and approaches zero only for asymptotically large T. In the vicinity of the deconfinement 
phase transformation, the dimensionless specific bulk viscosity (/s is expected to become 
large and even to diverge at a second-order phase transition 55H6JJ. However, there are 



examples of strongly coupled theories 62 



63], for which the contrary is found. 



ts for shear and bulk viscosities 

It 



Transport coefficients of strongly interacting systems have been calculated in a variety 
of different approaches. Despite the fact that lattice QCD studies of dynamical quantities, 
such as viscosities, are tedious to perform, they represent the first-principle non-perturbative 
method in the strong coupling regime. Recently, first resu 
obtained within lattice QCD simulations were reported in 

In quantum field theories, calculations of transport coefficients are often performed within 
the framework of linear response theory resulting in the known Kubo formulae 68M70| . This 
approach for determinin g sh ear and bulk viscosities was exploited for a weakly coupled 
scalar field theory in [tiI, 72]. Likewise, kinetic theory can be applied as a rigorous tool to 
quantify transport coefficients at weak coupling from a linearized Boltzmann equation. In 
fact, it was shown for weakly coupled scalar fields 71. 72 1 and for hot QED I73I, S] that 
both approaches yield equivalent results, when a linearized Boltzmann equation for excita- 
tions with temperature dependent self-energies and scattering amplitudes is considered. In 
QCD, kinetic-theory weak-coupling results were reported for the shear viscosity at leading 
logarithmic order 75] and at full leading order 76) in the QCD running coupling a s . These 
calculations found for a s < 0.25 a specific shear viscosity rj/s > 1. For the bulk viscosity 
perturbative QCD (pQCD) leading order results were presented in [77]. 
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Further investigations of transport coefficients utilize e. g. a specific ansatz for the spec- 
tral function entering the Kubo-expression for 77 ^7^8^ or apply sum rules for the spectral 
density of the trace of the energy-momentum tensor, thereby relating ( with basic ther- 
modynamic quantities 55H57I]. Moreover, within a parton cascade model based on kinetic 
theory 81] a value for rj/s ~ 0.5 was found [84], which is much smaller than the pQCD 
result rj/s ~ 2.7 [75|, |76| at comparable coupling a s ~ 0.1. Recent similar studies report 



rj/s ~ 0.8 |83[ and rj/s 



84j, 



851 ] instead. 



In this work, we study the properties of shear and bulk viscosities for the pure gluon 
plasma at finite temperature. The dynamics of the system is considered phenomenologically 
by assuming that it can be described in terms of gluonic quasiparticle excitations. In the un- 
derlying quasiparticle model (QPM), the (quasi)gluon excitations obey a medium-dependent 
dispersion relation, where the corresponding self-energy depends on T both explicitly but 
also implicitly via a phenomenological effective coupling G 2 (T) 86H89I]. In a pure gluon 
plasma, there are no conserved charges related to internal symmetries. Thus, T is the only 
thermal variable characterizing properties of the system. The QPM was successfully em- 
ployed to quantitatively describe lattice QCD results of equilibrium thermodynamics such 



as the EoS and related quantities [90 



94]. In the following, we apply this model to systems 



exhibiting small departures from thermal equilibrium by using the results of an effective 



71 



72( for shear and bulk viscosities derived for a quasiparticle 



kinetic theory approach 
system in [95]. 

We show that this quasiparticle picture combined with kinetic theory provides a powerful 
method to quantify the transport properties of the gluon plasma in a broad temperature 
range above the critical temperature for deconfinement T c . In particular, we show that in 
the high temperature limit the resulting transport coefficients exhibit parametrically the 
dependencies on T and the coupling known from pQCD. We also show that our results 
extrapolated into the non-perturbative regime are in agreement with available lattice QCD 
results. 

Our paper is organized as follows. In the next section, we recall the first-order non- 
equilibrium corrections to the energy-momentum tensor for the (quasi)gluon plasma and 
discuss features of its viscosity coefficients in relaxation time approximation as well as specify 
the employed quasiparticle parametrization. In sections IIHI and IIVt the shear and bulk 
viscosities are quantified and compared with available lattice QCD data. The summary of 
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our results is found in section [V] Throughout this work we use natural units, i. e. h = c — 
k B = I- 



II. SHEAR AND BULK VISCOSITIES 



For charge-neutral systems and in the absence of additional conservation laws, the local 
conservation equations for energy and momentum in relativistic hydrodynamics as effective 
theory of long-wavelength modes determine the evolution of the system. Phenomenologi- 
cally, the energy- momentum tensor T^, with equations of motion d v T^ v = 0, is defined as 
decomposition into a perfect fluid (or thermal equilibrium) part TffX and a part T^, which 



accounts for first-order dissipative effects [96l-l98j : T^ u = Tj^+Tff! with = eu^u" — PA 
and = QA^ u d a u a + rjS >1 ' / aj3 d 'u 15 , where e and P denote energy density and thermody- 
namic pressure of the system, respectively; both are related through the equation of state 
P = P(e). The projector A^ u = — u^u y with normalized four- velocity u^u^ = 1 and 



metric g^ u is orthogonal to the projector 



a. 3 



ACA" _|_ A^A V 



2A» v A af) /3 is also 



w^-orthogonal. The coefficients 77 and £ entering T%\ denote shear and bulk viscosities, 
respectivelv Note that if can be formulated solely as a function of e, then ( = is 



found 



99 



100|| . In the way T^ v is defined, one takes into account that from d^T^ = only 



four out of ten independent components of T^ u can be determined, while the others follow 
from expanding aro und in terms of small space-time gradients of the fluid four- 



velocity field m m 99|, llOOj . This expansion holds up to including further gradients or higher 



powers of u^. The Landau-Lifshitz condition [96 



97J e = T^UyUy = T^u^u u provides the 



only uniform definition of local flow for systems without additional conserved currents 44]. 



Basic definitions 



Concise expressions for the shear and bulk viscosities have been derived systematically 



111 



95[ within a quasiparticle framework. In relaxation time (r) approximation, one can 



deduce for the (quasi)gluon plasma 
1 f d 3 p 



V 



15T 



1 



— n(T)(l + crVT))-p 4 , 



C rj-i 



g n( T)(l + ^n(T))^{^-a]f-^ 



(1) 
(2) 



where n(T) = d{e E l T — with d — 16 is the equilibrium distribution function, dP/de 
denotes the squared speed of sound and a = T 2 dIi{T) / dT 2 . The excitations obey the disper- 
sion relation E 2 = p 2 + I1(T) with II(T) as temperature dependent self-energy. The equilib- 
rium distribution n(T) is a particular solution of an effective kinetic theory of Boltzmann- 
Vlasov type [96| with force F a {x) = {d a U(x))/ (2^/U{x)) |7l|, bj, satisfying p a F a = 
for the momentum p a ~ (E(p,x),p), in a state without entropy production where the 
collision term vanishes. Thermodynamic self-consistency is ensured, i.e. e(T) = T^u^u u , 
P(T) = -§T ( ^A Mi , and e = T 2 d(P(T)/T)/dT from e + P = Ts with s = OP/dT. 

Equ ation ([!]) for the shear viscosity, whic h was firs t derived in the quasiparticle framework 



in 



101[ | . coincides formally with the results 



102 



103| obtained within kinetic theory without 



the force term oc F a and without self-energy. Nonetheless, it differs from these results, as the 
effecti ye quasipartic le mass II (T) entering E reflects medium effects. This was also observed 



in 



60 



61 



104 



105]. 



The expression (T5]) for the bulk viscosity represents a special case of the general form of (, 



T I (0^)(1 + d~ l n{T))j^ {[E 2 - a] f - \p 2 } {X[E 2 



1^2 

3 



p 2 }, obtained by setting 
the momentum independent factor X = dP/de for invoking the Landau-Lifshitz condition. 

(Choosing instead * = (0n(T)(l + d^n^wl i - a ] f " W 2 }-^^ ~ a > 

would follow for (, which corrects a factor-2 mistake in the re sult deriv ed in 10l| .) Equa- 



tion (J2J) differs even formally from the previous results 1601 . 1611. Il02l - ll05| . This is because it 
contains the subtracted mass term (n(T) — a) as in 721. 1771 . 1 1 If ] and thus also the derivative 
of G 2 (T) with respect to T and not only II(T) alone. If one would use a temperature inde- 
pendent mass squared M 2 ins tead of II(T) in E, then both results for £ would correspond 



to the expression obtained in 



103 ]. 



Medium (i. e. EoS) effects influence the bulk viscosity more prominently than the shear 
viscosity. This is because £, just opposite to r], is relevant in processes, which change the 
volume of a system rather than its shape. Then, in the pure gluon plasma, the only way to 
relax locally the system to equilibrium is by changing the average energy of the excitations 
and, thus, the energy density. 

Discussing the general properties of ( from Eq. ([2]), we note that for a gas of particles 
with small and constant, i. e. temperature independent, mass a bulk viscosity ( ~ is found 
in the ultra-relativistic limit (p a p a — 0) because dP/de ~ 1/3. In the non-relativistic case, 



T, 



(i)i 



2T™ is obtained from the above decomposition of the energy-momentum tensor, 



6 



which in the local fluid rest frame vanishes due to the Landau-Lifshitz condition implying 
also C = 0- F° r massless particles one has exactly dP /de = 1/3, dll/dT = and E = \p\, 



98 



100 



102 



103 



106]. 



such that ( = is found from Eq. (J2]) in line with 

In the conformal limit, i. e. at asymptotically large T, the bulk viscosity vanishes, £ — > 0. 
This is because the squared factor in Eq. ([2]) can be rewritten as 

OP 1 _ 



[E 2 



p 

P la7 



1\ dP 

3 U^(n(T)-a) 



IdP ,dG^ 
2 de dT 2 



(3) 



for n(T) = T 2 G 2 (T)/2 (cf. section Hi Cj) . Consequently, in the combination (II (T) — a) the 
leading temperature dependence of IT(T) is canceled out, cf. also [771 ]. leaving solely the 
derivative of the effective coupling with respect to T. In the conformal limit, however, the 
only dimensionful scale is the temperature. This prohibits, in particular, a running of the 
coupling with T, i. e. dG 2 /dT = 0. Therefore, together with dP/de = 1/3 (cf. Eq. 
below), one finds 



P' 



dP 



IdP.dG 2 
2 97 #T2 



0. 



In other words, the terms (dP/de — 1/3) and (II(T) — a) are a measure for the deviation of 
the system from conformal invariance. 

B. Relaxation times 

To quantify shear and bulk viscosities from Eqs. Q and (J2]), respectively, one needs to 
specify r. Microscopically considered, the shear viscosity is dominantly determined by elas- 
tic two-body scatterings among excitations at typical momenta of order 0(T). The bulk 
viscosity, in contrast, measures relaxation of disturbances, which are caused by a slow uni- 
form expansion of the system. Consequently, it is determined by inelastic particle number 
changing processes and is sensitive to excitations with soft momenta 7l|, [72|, I75M77I. |98 |. 
This implies that the relaxation times for 7] and ( are in general different. In the following, 
however, we concentrate on elastic gluon-gluon 2 -h- 2 scatterings only. The correspond- 
ing relaxation rate, r _1 = h(va), is determined from the thermal-averaged total (elastic) 
scattering cross section (vcr) and the particle density in local equilibrium h, where v is the 
relative velocity between the two scattering particles. 
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The expression we employ as ansatz for r 1 originates from solving the Boltzmann equa- 



tion by a variational treatment including screening effects in the gluon-gluon scattering [107]. 
The result depends crucially on the ratio of maximum to minimum momentum transfer in 
the scattering process. The minimum momentum transfer is related to the Debye screening 
mass vhy). The maximum momentum transfer p max is instead limited by the typical particle 
momentum in the medium. For our quasiparticles this is 4 > p max /T > 3 in the temperature 
interval 1 < T/T c < 4, while for larger T, p m ax/T decreases extremely slowly with T and 
can be set constant. 

At lar ge te mperatures, the relaxation rate r _1 was found to be proportional to 
ln(l/m|,) 107]] . In order to extend this perturbative QCD result into the non-perturbative 
regime, we replace the running coupling g 2 = 47ra s entering tud by the effective coupling 
G 2 (T) from Eq. (jSJ) below. In this way, one arrives at the following phenomenological form 
for the inverse relaxation time 



T 



aiTG 4 ln(a 2 /G 2 



(4) 



The above expression is momentum independent and is assumed to be common for s hear and 



bulk viscosities. This may be considered as a rather crude approximation [75j, |77j, [108]. A 
more refined approach, employing a momentum dependent relaxation time including elastic 
and inelastic scatterings, was recently proposed in [95]. 



C. Temperature dependent self-energy 



For the quantitative discussion of Eqs. (ED) a nd (1 21), one also has to specify the self-energy 
n(T). We adopt the quasiparticle model 



effective coupling 



86 



-lasi 



G 2 {T) 



86H9JJ with II(T) = T 2 G 2 (T)/2 and use the 
16tt 2 



11 In 



T-T s 
T c /\ 



(5) 



inspired by the QCD 1-loop running coupling for purely gluonic matter. In fact, by con- 
struction, the effective coupling G 2 (T) reproduces g 2 (T) at large T, where a temperature 
dependence is introduced into g 2 by considering it at variable r enormalization point 2ttTc 



with c G [1/2,2], and as scale T c /A^- = 1.14 from 109l . IllOj is used. This implies, in 
particular, that G 2 (T) vanishes logarithmically as T — » oo. In Eq. (151), A and T s are two 
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T/T c 



FIG. 1: (Color online) Comparison of QPM results (solid curve) for the scaled interaction measure 
(e — 3P) / T as a function of T/T c for the pure gluon plasma with SU C (3) lattice QCD results 



from 



llll ] (boxes) and from 



112( | (circles). 



model parameters. 

This quasiparticle model works impressively well in describing lattice QCD results for the 
equation of state and related quantities by adjusting A and T s . Fig. [TJexhibits the comparison 
of QPM results for the scaled interaction measure (e — 3P)/T 4 as a functio n of T /T r with 



111 



the corresponding lattice QCD results obtained in pure SU C (3) gauge theory 
model parameters were chosen as A = 4.3 and T s = 0.73 T c with T c = 271 MeV 
order to obtain e and P [9l|, |92| , we have evaluated 

d 3 p 



1121. The 



113) . In 



e(T) 
P(T) 



(2vr) 3 
d 3 p p 2 



En(T)+B(T), 
n{T)-B{T), 



(6) 
(7) 



(2tt) 3 £ 3 

where B(T) is the solution of dB(T)/dU(T) — — | J j£M ' E n(T) , with integration constant 
B(T C ) = 0.19 T 6 4 . Thereby, it is implicitly assumed that B(T) depends on T only via n(T). 
As evident from Fig. [H both the maximum in the vicinity of T c , whic h sig nals the strong 
deviation of SU C (3) Yang-Mills theory from the conformal limit e = 3P jllj ]. as well as the 
decline towards the conformal limit for asymptotically large T observed in lattice QCD are 
well described within the QPM. The latter behaviour follows from the logarithmic decrease 
of the effective coupling G 2 (T) with T, while the position of the maximum is related to the 
inflection point of p/T A as a function of InT and depends on the choices made for A and T s . 



We note as an aside that the presence of d ynamical quarks significantly influences the 
thermodynamics of the system, cf. e. g. jl!5| for corresponding lattice QCD calculations 
using physical quark masses. Similarly, it is expected that also the transport coefficients are 
considerably different for QCD compared to the pure gluon plasma. It would be interesting 
to study how these differences eventually translate into the specific viscosity coefficients. 

The squared speed of sound as a function of T/T c following the QPM result for the 
interaction measure exhibited in Fig. [1] can be parametrized in the temperature interval 
T c <T <4T c by 



dP 



1 



tanh 



j2cn(T/T c -ir 



n=l 



(8) 



where the parameters read c\ = 4.78, C2 = —6.39, C3 = 4.60, C4 = —1.57 and C5 = 0.20. 
The quality of this fit is quantified by xVd.o.f. = 3.89 x 10~ 4 /16. The function dP/de is 
monotonously rising with T, from small values just above T c towards 1/3; at T = 4T C the 
squared speed of sound is only 2.6 % smaller than 1/3 . At T — > T c + , it decreases fastly to 
zero as expected for a first-order phase transition 116 ]. 

The high temperature expansion of P{T) from Eq. ([7]) for the pure SU C (3) plasma 



P(T) = aT A (1 - bG 2 (T) + C(G 4 )) 



(9) 



with a = d7T 2 /90 and b = 15/ (g?tt 2 ) was shown in 87j] to be in agreement with the pQCD 



result 



77| . The squared speed of sound at large T can be evaluated from standard thermo- 



dynamic relations dP/de = dP/de = (dP/dT)/(de/dT) = (dP/dT)/(Td 2 P/dT 2 ), where the 
individual derivatives read 



dP dC 2 
— = AaT 3 (l - bG 2 + 0(G 4 )) - abT 4 ^ + 0(dG A /dT) , 

rt 2 P AH 2 



(10) 



\2aT\\ - bG 2 + C(G 4 )) - 8abT 4 



0(d 2 G 2 /dT 2 ,dG 4 /dT). (11) 



dT2 — ■ -V" „ dT 

From Eq. ((5]) and at large T the following pattern is valid, 1 3> G 2 ^> T(dG 2 /dT) ^> 
T 2 {d 2 G 2 1 dT 2 ). Thus, the squared speed of sound in the gluon plasma at large temperatures 
is given by 

+ 0(G TdG /dT) . (12) 



de 3 36 dT 



The temperature derivative of the effective coupling follows from Eq. (jSJ) as 



(T-T s 



dG 2 
] 1T 



11 



57T- 



G\T) . 



(13) 



10 



Therefore, for large T > T s one finds T(dG 2 /dT) = -11G 4 (T)/(8tt 2 ) < 0. In other words, 
dP/de < 1/3 is obtained from Eq. ffl2|) . For asymptotically large T, dP/de approaches 1/3. 



III. THE SHEAR VISCOSITY 



In [76], the shear viscosity was obtained for the pure gfuon plasma with N c = 3 colours 
at next-to-leading log order of small running coupling g as 

ViT 3 

VNLL 9 nn[^/(gT)} 

with r]i = 27.126 and /x*/T = 2.765. With r from Eq. (J3J) and at large T, the expression for 
r\ in Eq. ([T]) exhibits parametrically the same dependencies on coupling and on temperature 
as the above perturbative result. Thus, within our formalism, by choosing the parameters 
in Eq. (j4]) appropriately, one can adjust the QPM to reproduce the above pQCD result well 
and then extrapolate the latter into the region close to T c . 

Taking into account the behaviour of G 2 (T) for large T, we choose a 2 = (/^*/T) 2 in order 
to recover the constant in the logarithm of ti^ll- Furthermore, the QPM expression for the 
shear viscosity from Eq. (TfJ can be written as rj = T A b(T) r. Here, b(T) is a dimensionless 
and with T monotonically rising function. It approaches slowly its asymptotic limit b(T — > 
00) = rfr(4)C(4)/(30vr 2 ) ~ 1.404 following the logarithmic decrease of G 2 (T) ->■ as T ->■ 00. 
Thus, by setting a\ ~ 2.587 ■ 1CT 2 , our r\ — > for T — > 00. Variations in a±, keeping 0,2 
fixed, allow to consider cases, where r/ — > t]nll for a large but finite temperature. 

The quantitative behaviour of the shear viscosity to entropy density ratio rj/s in the QPM 
as a function of T/T c is shown in Fig. |2]for different values of a\ and a 2 . The dashed curve 
corresponds to a\ ~ 2.587- 1CT 2 and a 2 = (/x*/T) 2 (Fit 1), while the solid curve corresponds 



to ai = 3.85- 10~ 2 and a 2 = 2(/^/T) 2 (Fit 2). The entropy density s in the QPM (cf. |9l|,|92j) 
was calculated by using the parametrization ()5]) for G 2 (T), entering also the viscosities, such 
that results on equilibrium thermodynamics for the pure gluon plasma obtained in lattice 
QCD are reproduced (see Fig. [TJ). The model results are compared with perturbative QCD 
results and with recent lattice QCD data in Fig. |5J 

As evident from Fig. |2l the ratio rj/s is very sensitive to the values of a\ and a 2 in 
Eq. (jlj). However, irrespective of the particular choices made for these parameters, the 
specific shear viscosity calculated from Eq. ([1]) decreases with decreasing T and exhibits 
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FIG. 2: (Color online) Specific shear viscosity rj/s as a function of T/T c . Dashed and solid curves 
exhibit the QPM results for Fit 1 and Fit 2, respectively (see text for details). The symbols denote 
lattice QCD data for pure S U r (3) gauge theory from jfjj] (full squares), from {(36] (diamonds 

open squares, open and full circles). The dotted line exhibits the 



and triangles) and from |6J 
Kovtun-Son-Starinets bound 



4514471] (t]/s)kss = l/(47r). The grey band at large T in dicates the 



117j]. The band 



perturbative QCD result t]nll/s with i]nll from [76] and entropy density s from 
was obtained by varying the renormalization point Ti in the QCD running coupling g(j£) in the 
range ttT < 7J < AirT. 

a minimum at T > T c . The precise location of this minimum is determined by the value 
of ct2- For even smaller T, rj/s increases again monotonically as T — > T c + . This picture 
is qualitatively in agreement with model predictio ns for the confined phase reporting a 



monotonically decreasing rj/s for T — > T c 118 



120]. For both fits, the ratio rj/s is found to 



be significantly smaller than unity in a temperature region around its minimum, which is in 
contrast to naive extrapolations of pure pQCD results. 

With the parameters from Fit 1, the specific shear viscosity approaches for asymptot- 
ically large T the result expected from perturbative calculations (see Fig. [2]). At smaller 



temperatures close 



QCD data [64 



66 



;o T c , agreement with some but not the general bulk of available lattice 
671 ] is found. Exhibiting a minimum of (rj/s) m i n ~ 0.6 at T ~ 1.75 T c , 



rj/s rises to quite large values for T — > T c + . This, h owever, s eems to be unfavoured by recent 



results obtained for T < T c in the confined phase 118Ml20 | 



Clearly, the extraction of transport coefficients from lattice QCD simulations is a dif- 
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T/T c 



FIG. 3: (Color online) Relaxation time r from Eq. for the parameters from Fit 1 (dashed curve) 
and Fit 2 (solid curve). See text for details. 



ficult task 
in 



121 



59 



79] and results seem far from being conclusive. Following the arguments 
122|, we focus on the lattice QCD results for r]/s from 641 ] . These results are 
best described with the parameters from Fit 2. In addition, as evident from Fig. [21 this 
parametrization shows fairly nic e ag reement with available lattice QCD results over a wide 
range of temperatures (see also [10 II ] ) . 

The behaviour of T]/s observed in Fig. [2] is a direct consequence of the particular form of 
r. In Fig. [3J the relaxation time is exhibited as a function of T/T c for the different parameter 
sets corresponding to Fit 1 and Fit 2. While being rather flat at moderate temperatures 
above T c , r increases rapidly for decreasing temperatures. In particular in the case of Fit 2, 
a pronounced rise of r in the vicinity of T c is observed. This behaviour is in line with the 
phenomenon of critical slowing down, indicating that the amount of time needed for relaxing 
energy and momentum density fluctuations increases rapidly near T c . 

Interestingly, the sharp increase of r close to T c is driven by the logarithm in Eq. ((4j). 
Without the logarithm, i. e. when assuming that r depends solely on (TG 4 ) -1 , the relaxation 
time would simply decrease with decreasing T in line with the increase of G 2 (T) for T — > T c + 
according to Eq. (JSJ). However, close to T c the factor In(a2/G 2 ) becomes dominant and 
determines the behaviour of r. The precise temperature dependence is cruci ally d riven by the 
parameter 02, which can lead to potentially huge quantitative differences 1081 ] . Increasing 



(decreasing) a2 results in a shift of the sharp increase in r to smaller (larger) T, apart from 



13 



an overall decrease (increase) of r. 

As evident from Eq. (j3J), the relaxation time diverges at T = T* for which G 2 (T*) = a%. 
For the parametrization of G 2 (T) employed in Fig. [TJ we find T* ~ 1.32 T c (Fit 1) and 
T* ~ 1.1 Tc (Fit 2). For T < T*, the relaxation time in Eq. becomes negative, indicating 
that the used ansatz is not valid anymore. 

At large T, the behaviour of r is dominated by its (1/T) dependence such that r — > 
for T —7- 00. However, in rj = T 4 6(T) r, this temperature dependence of r is compensated. 
The two different parametrizations Fit 1 and Fit 2 of r result in a finite difference between 
the corresponding shear viscosities denoted by 77 ^ for Fit 1 (using and ) and r]^ for 
Fit 2 (using and (4 ). At large T, the ratio rf^ /rf^ reads 



(i) 
(2) 



4! 



1 



In 2 



ln[(/i*/T) 2 /G 2 



,(1) 



1.49. 



Therefore, by construction, the ratio rf^Js approaches the pQCD result for asymptotically 
large T, whereas rj^/s does not reach this limit at any T, remaining at least a factor /a^ 
smaller. 

With the parameters from Fit 2, we find for the relaxation time r ~ 0.4fm/c for T > 
1.2 T c , whereas it increases sharply beyond 5fm/c as T — )■ T c + (see Fig. [3]). Similar values 



were used in the hydrodynamic simulations 



for discussing the influence of shear and 



bulk viscosities on observables measured in heavy-ion collisions such as the elliptic flow . 
Moreover, our results for r are numeric ally comp arable with results obtained in 78|, Il23| . 



where a different ansatz for r based on 



124 



125| was used reading 



2Nc 
8tt 



TG In 



2c 

7p 



(14) 



The parameter c is determin ed by dem anding r — > 00 as T — > T c , i. e. c = G 2 (T = T c )/2. 



The results for r]/s from 78|, Il04j . Il05l |. where this ansatz for r is employed, show a slightly 
flatter T dependence compared with our Fit 2. However, Eq. (1140 differs from Eq. (Hj) in 



its parametric dependence on G 2 . Therefore, as already discussed in 781 ]. the parametric 
dependence of rj on the coupling found in pQCD cannot be recovered with Eq. ffT4"j) for r -1 . 

As evident from Fig. HI the ratio r]/s for Fit 2 increases very mildly with increasing T. For 
T = 1.5T C , we find rj/s = 0.13, which is a factor 1/3 smaller than results obtained in 126 ] 
for the quark-gluon plasma with Nf = 3 degenerate quark flavours. At T = 3T C , rj/s ~ 0.29 
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T/T c 



FIG. 4: (Color online) Specific shear viscosity rj/s as a function of T/T c as in Fig. [2] but zoomed 
into the temperature region near T c . Symbols depict lattice QCD results as in Fig. [2] and the solid 
curve exhibits the QPM result for Fit 2. The area between the two grey dash-dotted curves shows 



the ex trapolatio n of the pQCD result t/tvll from 



QCD 



111 



761 ] scaled by the entropy density from lattice 



1 121 ] . where the renormalization point entering tjnll is varied between \i = nT and 



~p = 4vrT. 



0.096. A similar 



is still small and only a factor 3 larger than the minimum value (r)/s) m i 
value for the minimum was recently reported from a virial expansion approach 127 ] 



IV. THE BULK VISCOSITY 



The bulk viscosity exhibits different parametric dependencies on T and on the 



coupling 

in comparison with the shear viscosity as observed in pQCD calculations [75H77| . In the 
QPM at large T, one finds with the expression for the squared speed of sound Eq. (fl2"|) 
that Eq. ([3]), which enters ( quadratically, exhibits for thermal quasiparticle momenta the 
following structure 

dP 1 

~ 3 

Together with r _1 from Eq. 
\p\ ~ T, that 

c~J^g^. (16) 



V 



2 dt dT 2 



^T^p 2 - -T 3 ^ + 0(G Q ) ~ T 2 G\T) . (15) 
36 3T F 12 9T K ' K ' K ' 

this implies for the bulk viscosity at large T, where E ~ 
T 3 G 4 (T) 



ln(a 2 /G 2 (T)) 



15 



1 



10 

T/T c 



FIG. 5: (Color online) Specific bulk viscosity £/s as a function of T/T c . The so" 
the QPM result for Fit 2. It is compared with available lattice QCD data from 



id curve exhibits 
(full squares) 



65] 



and [6 jfl (open and full circles). The noticeable error bars for the results from 



651 ] represent 



recently reported from holographic QCD 



1281 ] 



conservative upper and lower bounds give n in [651 ]. The dotted curve shows for comparison results 



Keeping in mind the behaviour of G 2 (T) at large T, this result coincides with the leading log 
order parametric behaviour of £ obtained in pQCD calculations by an expansion in inverse 
powers of ln(l/a s ) jz3|. Nonetheless, our result for ( exhibits the parametric dependencies 
on T and G 2 (T) in Eq. (fT6l) only, because it contains the subtracted mass term (fl(T) — a) 
rather than the effective thermal mass fl(T) of the quasiparticle excitations alone. 

The QPM results for the specific bulk viscosity as a function of T/T c are depicted for 



Fit 2 in Fig. 0and compared with available lattice QCD results [65|, |67|. We find a positive 
but approximately vanishing ratio (/s for T > 2T C in line with lattice QCD. However, for 
T — > T 6 + , the specific bulk viscosity increases sharply following the behaviour of r. Similarly, 
for T — > T~ in the confine d phase, a monotonic increase of (/s was found in different model 



calculations 



118, 



129Ml3l[ ] . Both observations together suggest that for strongly interacting 



matter the ratio (/s develops a maximum in t he vic inity of T c 55M58I]. Such a behaviour is, 



132 ]. 



for instance, also seen in Lennard- Jones fluids 

Even though the QPM results are in qualitative agreement with lattice QCD, nonetheless, 
for T < 2T C they are systematically below the lattice QCD data 65j. However, the results 



from [65( are accompanied by noticeable uncertainties and a less rapid increase of (/s close 
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T/T r 



FIG. 6: (Color online) Combined specific sound channel (rj + 3£/4)/s as a function of T/T c . The 



QPM result for Fit 2 is sho wn b y the solid curve. It is compared with available 

i i rn 

taken either directly from 1221 ] (circles) or from combining results of [641 ] and 



attice QCD results 



651 ] (squares). The 



uncertainties in the latter arise from the upper and lower bounds in Q/ s depicted in Fig. [5j 



to T c seems to be likely 122]. In fact, refined calculations 1221 ] have been reported for 
the combined specific sound channel (77 + 3£/4)/s, however, only for temperatures, where 
the contribution from rj is dominant. A comparison of QPM predictions for the specific 
sound channel with available lattice QCD data from 1221 ] and from a combination of data 



from 



64 



65| is exhibited in Fig. |6l showing fairly nice agreement. 



Our results for (/s and r]/s exhibit both a pronounced behaviour close to T c , which might 
result in a sizeable influence of these transport coefficients on some observables measured 
in heavy-ion collisions. In fact, the observed possible rapid increase of rj/ s with de creasing 



T close to T c as well as the behaviour with T i n the con fined phase 



118 



refinements in viscous hydrodynamic simulations 



12fjj call for 



133H135] . The smallness of C/s over a 



wide range of T, in contrast, could suggest that the impact of bulk viscous effects is less 



important. Indeed, in |38j a small influence of (/s on the shear viscous suppression of the 
elliptic flow was found. Although Q/s becomes larger close to T c , it increases less than 
the corresponding relaxation time (see Figs. [3] and [5]) due to the thermodynamic integrals, 
cf. Eq. ([2]). In other words, critical slowing down hampers the influence of bulk viscosity 
as long as the corresponding (negative) bulk viscous pressure is small enough to ensure the 
stability conditions of viscous hydrodynamics [38]. This effect also tempers the tendency of 



17 



the fluid to becom e mechanically unstable against cavitation and clustering for large values 
ofC 



37 



133-138]. 



V. CONCLUSIONS 



Shear and bulk viscosities give important information about the transport properties 
of a medium. In this work, we calculated these transport coefficients for the pure gluon 
plasma within a quasiparticle model by assuming the plasma to be describable in terms of 
(quasi)gluon excitations with temperature dependent self-energy. 

The approach is based on an effective Boltzmann-Vlasov type kinetic theory for quasi- 
particle excitations with medium-dependent dispersion relation, which is consistent with 



the approach proposed in 



71 



721 ] . In local thermal equilibrium, this picture reduces to the 



quasiparticle model, which was employed successfully to describe lattice QCD equilibrium 



thermodynamics in [90|, l92H94j . Thermodynamic self-consistency in this model is a direct 



consequence of the consistency condition that must be imposed in order to fulfill energy and 
momentum conservation in the kinetic theory 



95 



101J. We used expressions for shear and 



bulk viscosities derived within relaxation time approximation to the kinetic theory. In the 
shear viscosity, medium effects appear only implicitly via the dispersion relation. In the bulk 
viscosity, medium effects are in addition explicitly reflected by the subtracted mass term, 
which contains the derivative of the effective coupling with respect to the temperature. 

Including only contributions from elastic gluon-gluon scatterings to the relaxation time, 
we have shown that at large T our results for the transport coefficients exhibit the same 
parametric dependencies on temperature and coupling as found in pQCD calculations 



75 



771 ] . For this to hold true, it is crucial that the bulk viscosity contains the subtracted mass 
term rather than the self-energy alone. 

In the temperature region close to T c , fairly nice quantitative agreement with available 
lattice QCD results for the transport coefficients is found. The specific shear viscosity 77/s 
exhibits a minimum in the vicinity of T c , while the specific bulk viscosity (/s, being large 
close to T c , falls off rapidly to its conformal limit for larger T. 

For the specific shear viscosity, we find a rather mild increase with temperature. In 
fact, the ratio r\j s is smaller than unity even at temperatures being a few times larger than 
T c , e. g. at T ~ 3T C it is still about 0.3. Even though our investigations were limited to 



18 



the case of a pure gluon plasma, one might expect that for thermal conditions relevant in 
heavy- ion collisions at the LHC (Large Hadron Collider), the ratio rj/s is small such that 
iC h eavy-ion programme based on ideal hydrodynamic simulations are 
14o| . However, to make definite predictions it would be of importance 
to include also quark degrees of freedom into our quasiparticle approach. 



predictions tor t he L 



meaningful 192 



139 
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